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The energy and momentum densities of the fields of a free 
electron in a plane electromagnetic wave include interference 
terms that are the classical version of the "dressing" of the 
electron the arises in a quantum analysis. The transverse me- 
chanical momentum of the oscillating electron is balanced by 
the field momentum resulting from the interference between 
the driving wave and the static part of the electron's field. 
The interference between the wave and the oscillating part 
of the electron's field leads to a longitudinal field momentum 
and a negative field energy that compensate for the longi- 
tudinal momentum and kinetic energy of the electron. The 
interference terms are dominated by the near zone, so that as 
the wave passes the electron by the latter reverts to its energy 
and momentum prior to the arrival of the wave. 



I. INTRODUCTION 

The behavior of a free electron in a electromagnetic 
wave is one of the most commonly discussed topics in 
classical electromagnetism. Yet, several basic issues re- 
main to be clarified. These relate to the question: to 
what extent can net energy be transferred from an elec- 
tromagnetic pulse (such as that of a laser) in vacuum to 
a free electron? 

These issues are made more complex by quantum con- 
siderations, including the role of the "quasimomentum" 
of an electron that is "dressed" by an electromagnetic 
wave 0. 

As a small step towards understanding of the larger is- 
sues, we consider a simpler question here. The response 
of a free electron to a plane electromagnetic wave is oscil- 
latory motion in the plane perpendicular to the direction 
of the wave, in the first approximation. Thus, the elec- 
tron has momentum transverse to the direction of the 
wave. However, the wave contains momentum only in 
its direction, and the radiated wave contains no net mo- 
mentum (in the nonrelativistic limit). How is momentum 
conserved in this process? 

The general sense of the answer has been given by 
Poynting Q , who noted that an electromagnetic field can 
be said to contain a flux of energy (energy per unit area 
per unit time) given by 



magnetic field (taken to be in vacuum throughout this 
paper) and c is the speed of light. 

Poincare |j) noted that this flow of energy can also be 
associated with a momentum density given by 
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Hence, in the problem of a free electron in a plane electro- 
magnetic wave we are led to seek an electromagnetic field 
momentum that is equal and opposite to the mechanical 
momentum of the electron. 

In this paper we demonstrate that indeed the mechan- 
ical momentum of the oscillating electron is balanced by 
the field momentum in the interference term between the 
incident wave and the static field of the electron. We are 
left with some subtleties when we consider the interfer- 
ence between the incident wave and the oscillating field 
of the electron. 



II. GENERALITIES 

A. Motion of an Electron in a Plane Wave 

We consider a plane electromagnetic wave that prop- 
agates in the +z direction of a rectangular coordinate 
system. A fairly general form of this wave is 



xEjc cos(fcz — tot) — yE y sin(fcz — cut), 
<J <j ' 



B waV c = ~x-E y sin(fcz — cut) + yE x cos(kz — cut), (3) 
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where to — kc is the angular frequency of the wave, k = 
27r/A is the wave number and x is a unit vector in the 
x direction, etc. When either E x or E y is zero we have 
a linearly polarized wave, while for E x = ±E y we have 
circular polarization. 

A free electron of mass m oscillates in this field such 
that its average position is at the origin. This simple 
statement hides the subtlety that our frame of reference 
is not the lab frame of an electron that is initially at rest 
but which is overtaken by a wave HIJ. If the velocity of 
the oscillating electron is small, we can ignore the v/cx B 
force and take the motion to be entirely in the plane z = 
0. Then, (also ignoring radiation damping) the equation 
of motion of the electron is 



in Gaussian units, where E is the electric field, B is the 



mx = eE wavc (0, t) = e(xE x coscot + yE y sincji). (4) 
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Using eq. (||) we find the position of the electron to be 



■ (ptEx cos cot + yE v sin Lot) , 



(5) 



and the mechanical momentum of the electron is 



Pmech = ItX = —(?C.E X SUlLOt — yE y COS Lot) . (6) 



The root-mean-square (rms) velocity of the electron is 



Wrms = V (x 2 + V 2 ) = \ 

muj V 
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The condition that the v/c x B force be small is then 
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where the dimensionless measure of field strength, r/, is 
a Lorentz invariant. Similarly, the rms departure of the 
electron from the origin is 



eE r , 
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Thus, condition (|J) also insures that the extent of the mo- 
tion of the electron is small compared to a wavelength, 
and so we may use the dipole approximation when con- 
sidering the fields of the oscillating electron. 

In the weak-field approximation, we can now use (|^) 
for the velocity to evaluate the second term of the Lorentz 
force: 
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This term vanishes for circular polarization, in which case 
the motion is wholcly in the transverse plane. However, 
for linear polarization the v/c x B force leads to oscilla- 
tions along the z axis at frequency 2lo, as first analyzed 
in general by Landau 0. For polarization along the x 
axis, the x-z motion has the form of a "figure 8" , which 
for weak fields (77 <C 1) is described by 



eE x 
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If the electron had been at rest before the arrival of 
the plane wave, then inside the wave it would move with 
an average drift velocity given by 



V 2 /2 
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along the direction of the wave vector, as first deduced 
by McMillan ||. In the present paper we work in the 
frame in which the electron has no average velocity along 
the z axis. Therefore, prior to its encounter with the 
plane wave the electron had been moving in the negative 
z direction with speed given by (|l2|). 



B. Field Momentum 

The fields associated with the electron can be regarded 
as the superposition of those of an electron at rest at 
the origin plus those of a dipole consisting of the actual 
oscillating electron and a positron at rest at the origin. 
Thus, we can write the electric field of the electron as 
Estatic + E osc and the magnetic field as B osc , where the 
oscillating fields have the pure frequency lo in the low- 
velocity limit. 

The entire electromagnetic momentum density can 
then be written 
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However, in seeking the field momentum that opposes 
the mechanical momentum of the electron, we should 
not include either of the self-momenta E wavo x B wavo or 
(E stat ic + E osc ) x B osc . The former is independent of the 
electron, while the latter can be considered as a part of 
the mechanical momentum of the electron according to 
the concept of "renormalization" . 

We therefore restrict our attention to the interaction 
field momentum 
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We recall from eqs. (|6|) and (gjj) that the transverseme- 
chanical momentum of the oscillating electron has pure 
frequency lo. Since the wave and the oscillating part 
of the electron's field each have frequency lo, the term 
contains harmonic functions of lo 2 , which can 



wave, osc 



be resolved into a static term plus ones in frequency 
2lo. Hence we should not expect this term to cancel 
the mechanical momentum. Rather, we look to the term 
Pwavc. static, since this has pure frequency lo. 



III. THE MOMENTUM tr WAVE, STATIC 

The static field of the electron at the origin is, in rect- 
angular coordinates, 



E 



static 
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where r is the distance from the origin to the point of 
observation. Combing this with eq. (||) we have 
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+yzE y sin(fcz — Lot) (18) 
+z[xE x cos(kz — tot) — yE y cos(fcz — wt)]}. 

When we integrate this over all space to find the total 
field momentum, the terms in z vanish as they are odd 
in either x or y. Likewise, after expanding the cosine and 
sine of kz — tot, the terms proportional to z cos kz vanish 
on integration. The remaining terms are thus 



(19) 
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after an elementary volume integration. 

It is noteworthy that the integration is independent 
of any hypothesis as to the size of a classical elec- 
tron. Indeed, the integrand of ([l9|) can be expressed 
as cos0sin(fcrcos6')/r 2 via the substitution z = rcos9. 
Hence, the integral over a spherical shell is independent 
of r for kr <C 1 , and significant contributions to the inte- 
gral occur for radii up to one wavelength of the electro- 
magnetic wave. This contrasts with the self-momentum 
density of the electron which is formally divergent; if the 
integration is cut off at a minimum radius (the classical 
electron radius), the dominant contribution occurs within 
twice that radius. 

Thus, we have demonstrated the principal result of this 
paper. 



IV. THE MOMENTUM rwAVE,OSC 

Several subtleties in the argument appear when we 
consider the other interference term in the momentum 
density ([bi]). For this we must first display the electro- 
magnetic fields of an oscillating electron. 



A. The Fields E osc and B osc 

Since we restrict our attention to an electron that os- 
cillates with amplitude much less than a wavelength of 
the driving wave, and the electron attains velocities that 
are much less than the speed of light, it is sufficient to 
use the dipole approximation to the fields of the electron. 
While these fields are well known, they are typically pre- 
sented in imaginary notation, of which only the real part 
has physical significance. This notation is very useful 
for discussions in which only time-averaged behavior is 
of interest. However, we wish to consider the details of 
momentum balance at an arbitrary moment, and it is 
preferable to use purely real notation. 



We begin by noting that the retarded vector potential 
of the oscillating electron at a point r at time t can be 
written 



-A. sc(r, ^) 



e ±{t' = t - r/c) 
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[xE x sin(fcr — tot) + yE y cos(fcr — tot)], 



using eq. (||) for the motion x of the electron. The oscil- 
lating part of the scalar potential is obtained by integra- 
tion of the Lorentz gauge condition: 
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We find 
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The constant static potential is omitted in the above. 

The scalar potential could also be deduced from the 
retarded potential of a moving charge. Equation ( |2^ ) 
results on expanding the retarded distance to first order 
in the field strength of the plane wave. 

The electric and magnetic fields are, of course, found 
from the potentials via 



B = V x A 



and 
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The lengthy expressions for the rectangular components 
of the fields are 
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These expressions can also be deduced from the Lienard- 
Wiechert forms for the fields of an accelerated charge, 
keeping terms only to first order in the strength of the 
plane wave. 



B. Components of P 
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Since the wave fields have no z component, the x com- 
ponent of P wa ve,osc is given by 
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From eqs. ([24|) and ( |25| ) we see that both B osc ^ z and E osc ^ z 
are odd in either x or y. Therefore, the volume integral 
of -P W avc,osc,K vanishes, and we do not consider it fur- 
ther. Likewise, P W nvc.osc.y vanishes on integration. This 
confirms the claim made at the end of sec. II that the 
interference term P wa vo,osc is not relevant to the balance 
of transverse momentum between the electron and the 
fields. 

However, the z component of P W avo,osc does not vanish 
on integration, and requires further discussion. As the 
details include some surprises (to the author) I present 
them at length. 

-*wave,osc,z = 
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The terms of -P W avc,osc,z that are proportional to E y E y 
are odd on both x and y, and so will vanish on integra- 
tion. 

We now consider the implications of eq. (^7j) separately 
for waves of circular and linear polarization. 



C. Circular Polarization 

For a circularly polarized wave, we have E 2 — E 2 . 
Consequently the dimensionless measure of field strength 
is r\ = eE x /mLoc — eEy/mujc, according to (||). The 
prefactors e 2 E 2 / Airmui 2 c and e 2 E 2 / Airmui 2 c can there- 
fore both be written rfmc/Ai:, and have dimensions of 
momentum. 

The terms of eq. (|27]) in E 2 and E 2 can be combined 
in pairs via the identities 

cos(fcz — Lot) cos(fcr — Lot) + sin(fcz — Lot) sin(fcr — Lot) 
= cos kz cos kr + sin kz sin kr, (28) 

and 

sin(/cz — cot) cos(kr — Lot) — cos(fcz — Lot) sin(fcr — Lot) 
— sin kz cos kr — cos kz sin kr. (29) 

A detail: the second term of eq. (|27| ) in E 2 contains fac- 
tors of x 2 , while second term of in E 2 contains factors 
of y 2 . But during integration, we can replace y 2 by x 2 , 
after which the terms can be combined via ( 28 - 29|) . 

We see already that the volume integral of -P W ave,osc,z 
will contain no time dependence! 

On integration, terms such as f(x,r)sinkz and 
g(x, r)z cos kz that are odd in z will vanish. The inte- 
grated field momentum is thus, 
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where 1\ is the volume integral whose integrand is 



4 



k 2 z 



sin kz sinfcr - 



kz 



■ sin kz cos kr 



*$kx k . 

— - A ) cos kz sin kr 



r 



fc 2 y 2 



— = y cos kz cos fer. 



(31) 



We return to the significance of eq. (|30|) after describ- 
ing the evaluation of integral I\ . 

As seen from eq. (|30|), the integral I\ must be dimen- 
sionless, although it is apparently a function of the wave 
number k. However, the form of ([H]) indicates that I\ 
is actually independent of the length scale, so we can set 
k = 1 during integration. 

To perform the integration we consider a volume el- 
ement r 2 dr dcosO dip in a spherical coordinate system 
with angle 9 defined relative to the z axis. It is more 
convenient to keep z = r cos 9 as a variable of integra- 
tion, using dz — rd cos 9. Then the volume integration 
has the form 



1=1 rdr / dz I 

Jv JO J-r JO 



(32) 



Most terms of (31) are independent of so their 
integral is just 2ir. For the terms in x 2 , we have 



2 • 2 

r sin i 



cos 



ir{r 2 - z 2 ). (33) 



While each of the four main terms of (|3l|) diverges 
on integration, it turns out that the two terms in cos z 
taken together are finite (and likewise for the two terms 
in sinz). We find that 
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where 
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From detailed evaluation of the radial integral, we find 
that the integrand approaches a constant value as r goes 
to zero, and that the contribution to the integral at large 
r diminishes as 1/r. That is, the principal contribution 
is from the region kr w 1. 

We are left with the result (^o| ) that the integral of the 
interference term in the field momentum density has a 
constant longitudinal term for an electron oscillating in 
a circularly polarized wave. 

Recall that we have performed the analysis in a frame 
in which the electron has no longitudinal momentum. 
However, as remarked in sec. IIA, prior to its encounter 
with the wave, the electron had velocity v z = —i] 2 c/2 (as- 
suming rj 2 <C 1), and therefore had initial mechanical mo- 
mentum p m ech,z = ^rj 2 mc/2. So, we would expect that 
this initial mechanical momentum had been converted to 
field momentum, if momentum is to be conserved. 

The result (]3(j) can be described as a kind of "hidden 
momentum" fflL whose appearance can be surprising if 
one ignores the physical processes needed to arrive at the 
nominal conditions of the problem. 

We continue to be puzzled as to why the result (|3(]) is 
8/3 times larger than that required to satisfy momentum 
conservation. 



D. Linear Polarization 

Consider now the case of a linearly polarized wave with 
electric field along the x axis. Then E lms = E x /y/2, and 
the prefactors in ( p7| ) can be written as r) 2 mc/2ir. 

The remaining terms in the momentum density 
-Pwave,osc.z have time dependences that can be expressed 
as sums of pure frequencies via the identities 



2 cos(fcz — Lot) cos(fcr — cot) 
— cos kz cos kr + sin kz sin kr 

+(cos kz cos kr — sin kz sin kr) cos 2tot 
+(cos kz sin kr + sin kz sin kr) sin 2ujt, 
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and 



2 cos(fcz — tot) sin(fcr — cut) 
cos kz sin kr — sin kz cos kr 

+(cos kz sin kr + sin kz cos kr) cos 2tot (38) 
+(sin kz sin kr — cos kz cos kr) sin 2u>t, 



Inserting these into eq. ( |27| ) and keeping only those 
terms that are even in z, we find the integrated field 
momentum to be 



Pwave,osc,2i — / -fwavc,osc,2i 
IV 



rj 2 mc 
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(Ji + I 2 cos 2tot + I 3 sin 2u)t), (39) 



(36) 

where integral 1\ = 167r/3 has been discussed in (||l]-[56|), 
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h = -Ia + Ib =~y> 
and integral -Z3 has the integrand, 
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On evaluation, ^3 = 0. 

Hence, the longitudinal component of the interference 
field momentum of a free electron in a linearly polarized 
wave is 



P 



wave . osc , z 



1 2 

-n mc - 



-r\mc cos 2u>t. 
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The constant term is the same as that found in eq. < p(\ ) 
for circular polarization, and represents the initial me- 
chanical momentum of the electron that became stored 
in the electromagnetic field once the electron became im- 
mersed in the wave. 

As for the second term of (f42|), recall from eq. pi] ) 
that for linear polarization the electron oscillates along 
the z axis at frequency 2lo. Hence the z component of 
the mechanical momentum of the electron is 



Pmcch.; 



r) 2 mc 



cos 2tot. 



(43) 



The term in p W avc,osc,z at frequency 2lo is —4/3 of the lon- 
gitudinal component of the mechanical momentum asso- 
ciated with the "figure 8" motion of the electron. Thus, 
we have not been completely successful in accounting for 
momentum conservation when the small, oscillatory lon- 
gitudinal momentum is considered. 

The factors of 4/3 and 8/3 are presumably not the 
same as the famous factor of 4/3 that arise in analyses of 
the electromagnetic energy and momentum of the self 
fields of an electron jl0|,[Ll]]. A further appearance of 
a factor of 8/3 in the present example occurs when we 
consider the field energy of the interference terms. 



V. THE INTERFERENCE FIELD ENERGY 

It is also interesting to examine the electromagnetic 
field energy of an electron in a plane wave. As for the 
momentum density (13), we can write 
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for the field energy density. Again, we no not consider 
the divergent energies of the self fields, but only the in- 
terference terms, 
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In general, the interference field energy density is os- 
cillating. Here, we look for terms that are nonzero after 
averaging over time. We see at once that 
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(48) 



since all terms have time dependence of cos tot or sinwt. 
In contrast, ([/wave, osc) will be nonzero as its terms are 
products of sines and cosines: 
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The terms in E x E y will vanish on integration over vol- 
ume. The various time averages are 

(2 cos(fcz — Lot) cos(fcr — tot)) 

= cos kz cos kr + sin kz sin kr. 
(2 sin(fcz — tot) cos(fcr — tot)) 

= sin kz cos kr — cos kz sin kr, 
(2 cos(fcz — tot) sin(fcr — Lot)) 

— cos kz sin kr — sin kz cos kr, 
(2 sin(fcz — Lot) sin(fcr — Lot)) 

= cos kz cos kr + sin kz sin kr. (50) 
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After performing the time average on eq. (J49|), we keep 
only terms that are even in z. These terms have the form 
(pl|), and so we find that 



Mint = / (Uv. 
Jv 



e 2 {El + El 
$>innuj 2 



T 2. 2 

-ll = rj mc 



(51) 

for waves of either linear or circular polarization. As 
with the case of the interference field momentum, this 
interference field energy is distributed over a volume of 
order a cubic wavelength around the electron. Being an 
interference term, its sign can be negative. 
We can interpret the quantity, 
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as compensation for the relativistic mass increase of the 
oscillating electron, which scales as d^/c 2 and hence as 
rj 2 (for small rj, recall eq. (0)). Indeed, a general result 
for the motion of an electron in a plane wave of arbitrary 
strength rj is that its rms relativistic mass, often called 
its effective mass, is 



Weg = my 1 + r\ 2 . 
For small rj, the increase in mass is 

Am 



1 2 
—7) m. 

2 ' 



(53) 



(54) 



Thus, the decrease in field energy due to the interfer- 
ence terms between the electromagnetic fields of the wave 
and electron is —8/3 times the mass increase it should 
compensate. 



VI. DISCUSSION 

A. Temporary Acceleration 

We remarked in sec. IIA that the preceding analysis 
holds in the average rest frame of the electron. If instead 
the electron had been at rest prior to the arrival of the 
plane wave, the velocity of the average rest frame would 
be v z = (t7 2 /2)/(1 + 77 2 /2). For this, the amplitude of the 
plane wave is presumed to have a slow rise from zero to 
a long plateau at strength tj, followed by a slow decline 
back to zero. 

Once the wave has passed by the electron, the inter- 
ference field energy, (|l]), goes to zero since the integral 
is dominated by the contribution at distances of order a 
wavelength from the electron. Hence, the electron's ki- 
netic energy must return to zero (or to its initial value 
if that was nonzero). A plane wave, or more precisely, 
a long pulse that is very nearly a plane wave, cannot 
transfer net energy to an electron. The acceleration of 
the electron from zero velocity to v z is only temporary, 
i.e., for the duration of the plane wave pulse. 

This result was first deduced by di Francia [Q and by 
Kibble by different arguments. 



B. The Radiation Reaction 

Our analysis of the energy balance of an electron in 
a plane wave is not quite complete. We have neglected 
the energy radiated by the electron. Since the rate of 
radiation is constant (once the electron is inside the plane 
wave), the total radiated energy grows linearly with time, 
and eventually becomes large. The interference energy, 
(pl|), is constant in time, and hence cannot account for 
the radiated energy. 

More to follow 



VII. APPENDIX: LIENARD-WIECHERT FIELDS 

As an alternative to the dipole approximation, we con- 
sider the use of the Lienard-Wiechert potentials and fields 
of a moving electron. We have limited our analysis to the 
case of a weak plane wave (rj <Cl), for which the velocity 
of the electron is always small (/3 = «/c<Cl). In this 
case we may approximate the time-dependent part of the 
fields of the electron as proportional to the strength of 
the field of the plane wave (proportional to rj. Then we 
find that the Lienard-Wiechert fields of the electron are 
the same as the fields in the dipole approximation. 

We can show this in two ways. First, we verify that the 
Lienard-Wiechert potentials reduce to eqs. (|2^) and (|2^). 
Second, we can verify directly that the Lienard-Wiechert 
fields are the same as eqs. (p3) and (25). 



The Lienard-Wiechert potentials are 



R(l-P-h 



A = 



e/3 



i?(l-/?-n 



(55) 



where the electron is at postion x, the observer is at r, 
their separation is R = r x, the unit vector n is R/i?, 
and the brackets, [ ], indicate that quantities within are 
to be evaluated at the retarded time, t' = t — R/c. 

We work in the average rest frame of the electron. In 
the weak-field approximation we ignore the longitudinal 
motion of the electron, (|ll|), which is quadratic in the 
strength of the plane wave. Then the velocity vector of 
the electron is 



(3(t) = (x-E x smut — yE y cos tot) . 



(56) 



from eq. (|6|). The retarded velocity is thus, 
[(3] = P{t' = t - R/c) 



(57) 



(xE x sm(kR — ujt) + yE y cos(fcr — uit)) 



mujc 

Distance R differs from r because the electron's oscilla- 
tory motion takes it away from the origin. However, the 
amplitude of the motion is proportional to strength of the 
plane wave. Hence, we may replace R by r in eq. ( [5^ ) 
with error only in the second order of field strength. 

Since the vector potential includes a factor (3 in the 
numerator, we can replace R by r and 1 — (3 ■ n by 1 
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in the first order in the field strength of the plane wave. 
Thus, 



A 



rriLocr 



(5tE x sin(fcr — tot) + yE y cos(fcr — tut)) . 



(58) 



in agreement with eq. (p3|). 

In the scaler potential, we first bring [3 to the numer- 
ator: 



e[l + ■ n] 
[R] ' 



Unit vector [n] differs from unit vector f due to the oscil- 
lation of the electron, which is proportional to the field 
strength of the plane wave. For the scalar potential, how- 
ever, we must expand the factor 1/[R] to first order in 
the field strength. Now, 



[R] = |r - x(t')| = v/r 2 -2r-x(t')+x 2 (t'), (60) 



with 



x(i') « T (iE x cos Lot' + E v cos cot') (61) 

e 



moo* 



(xE x cos Lo(kr — Lot) — yE y cos(fcr — cut)) , 



again approximating R by r in the arguments of the co- 
sine and sine, accurate to first order in the field strength. 
Hence, 



1^(1 + -<«'» 



1 f (xE x cos(fcr — Lot) — yEy sin(/cr — uit)) 



r 



Altogether, 



r muj 

TP < k y 



E. r 



kx 



sin(fcr — Lot) H — - cos(fcr — Lot) 

!J*£ Jit) 



y 

cos(fcr — Lot) sva{kr — Lot) 
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(62) 



(63) 



in agreement with eq. (|22|). 

Similarly, we could proceed from the Lienard-Wiechert 
fields, 



E = e 



n-/3 



_ 7 2 (l-/3-n) 3 i? 2 
B = [n x E]. 



nx Un-P) x j3 1 
(l-[3-h) 3 R 



(64) 



After some work, we find that these fields are the same 
as eqs. (|24]-p5|), to first order in the strength of the plane 
wave. 
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